Phase shift induces currents in a periodic tube 
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Abstract 

The average current of an overdamped Brownian particle moving along the axis of a three- 
dimensional periodic tube is investigated in the presence of a symmetric potential and a temporally 
symmetric unbiased external force. Reduction of the spatial dimensionality from two or three 
physical dimensions to an effective one dimensional system entails the appearance of not only an 
entropic barrier but also an effective diffusion coefficient. We find that the phase shift between 
the tube shape and the potential can break the symmetry of the effective potential and can induce 
net currents. Under optimal condition, the current as a function of temperature and the phase 
shift possesses many extrema of alternating signs. The current may reverse its direction several 
times when temperature or the phase shift is changed. Our model is to describe the movement of 
molecular motors along microtubule. 
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I. INTRODUCTION 

Molecular motors are protein molecules that can convert chemical energy, usually in the 
form of adenosine triphosphate (ATP), into mechanical forces and motion. Most organisms 
have many different motors that are specialized for particular purposes such as cell division, 
cell crawling, cell shape maintenance and movements of internal organellesi*^. Three dif- 
ferent linear families of molecular motors have been identified^: Kinesins and dyneins that 
move along microtubule, myosins that move along actin filaments. From a theoretical point 
of view, molecular motors are microscopic objects that move unidirectionally along periodic 
structures. The problem of this unidirectionality belongs to a larger class of such problems 
involving rectifying processes at small scale. Noise-induced transport phenomena play a 
crucial role in molecular motors^"-^. A molecular motor is usually described by a ratchet 
system which is generally defined as a system that is able to transport particles in a periodic 
structure with nonzero macroscopic velocity in the absence of macroscopic force on average. 

In these systems, directed Brownian motion of particles is generated by nonequilibrium 
noises in the absence of any net macroscopic forces and potential gradients. Typical ex- 
amples are rocking ratchets--, hashing ratchets-, diffusion ratchets^fi, correlation ratchets^i^ 
and white-shot-noise ratchets^. In all these studies, the potential is asymmetric in space. It 
has shown that a unidirectional current can also appear for spatially symmetric potentials if 
the external random force is either asymmetric or spatially dependent. If spatially periodic 
structures are exposed to an additive Poissonian white shot noise, a macroscopic current 
occurs even in the absence of spatial asymmetry^. 

Most studies have revolved around the energy barrier. The nature of the barrier depends 
on which thermodynamic potential (internal energy or Helmholtz free energy) varies when 
passing from one well to the other, and its presence plays an important role in the dynamics 
of the system. Whereas energy barriers are more frequent in problem of solid-state physics 
(metals and semiconductors, coupled Josephson junction and photon crystal), entropy barri- 
ers are often encountered in soft-condensed matter and biological systems. Entropy barriers 
may appear when coarsening the description of a complex system in order to simplify its 
dynamics. Reguera and co-workers^^ used the mesoscopic nonequilibrium thermodynamics 
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theory to derive the general kinetic equation of the motor system and analyzed in detail the 
case of diffusion in a domain of irregular geometry in which the presence of the boundaries 
induces an entropy barrier when approaching the dynamics by a coarsening of the descrip- 
tion. In their recent work^^, they studied the current and the diffusion of a Brownian particle 
moving in a symmetric channel with a biased external force. They found that temperature 
dictates the strength of the entropic potential, and thus an increasing of temperature leads 
to a reduction of the current. 

The focus of previous works on current were limited to the case of energy barriers in the 
presence of asymmetry in space or time. The present work is extended to the study to case 
of entropic barriers in the presence symmetry in space and time. We emphasize on finding 
how phase shift between the tube shape and the potential can induce a net current. Our 
model is to describe the movement of kinesins and dyneins along microtubule. 

II. THE CURRENT INDUCED BY A PHASE SHIFT 

Phenomena of current reversals may be of interest in molecular motors^^. kinesins and 
dyneins move along tubulin filaments towards their plus and minus extremities, respectively. 
It is well known that the two current reversals effect allows one pair of motor proteins (two 
members)to move simultaneously in opposite directions along the microtubule inside the 
eukaryotic cells. One member of the pair moves one way along the tube while the other 
member of the pair moves in the opposite direction. In order to describe the movement 
of kinesins and dyneins along tubulin, a periodic tube model is proposed. We consider 
a Brownian particle moving in a symmetric periodic tube [Fig. 1] in the presence of a 
temporally symmetric unbiased external force. Its overdamped dynamics is described by 
the following Langevin equations written in a dimensionless form^^*^. 



where x, z are the three-dimensional (3D) coordinates, rj is the friction coefficient of 
the particle, is the Boltzmann constant, T is the absolute temperature and C,x,y,z{t) 




(1) 




(2) 




(3) 
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is the Gaussian white noise with zero mean and correlation function: < ^i{t)^j{t') >= 
26ij5{t — t') for i,j = x,y,z. < ... > denotes an ensemble average over the distribution 
of noise. 6{t) is the Dirac delta function. Imposing reflecting boundary conditions in the 
transverse direction ensures the conflnement of the dynamics within the tube, while periodic 
boundary conditions are enforced along the longitudinal direction for the reasons noted 
above. F{t) is a temporally symmetric unbiased external force along the x direction which 
is imparted as the result of energy gained via ATP hydrolysis and satisfles^""- 

[ -Fo, nr + ^T <t<{n + l)r, 

where r is the period of the unbiased force and Fq is its magnitude. 

U{x) is a symmetric potential with periodicity L which is from the track of molecular 
motors [Fig. 1], 

U(x)=Qsm{'^), (5) 

where Q is the amplitude of the potential. The shape of the tube (describing the structure 
of molecular motors) is described by its radius, 

Co'(x) = a sin(— — \- (p) + b, (6) 
ij 

where a is the parameter that controls the slope of the tube, is the relative phase shift 
between u){x) and U{x). The radius at the bottleneck is b — a. 

The movement equation of a Brownian particle moving along the axis of the 3D (or 
2D) tube can be described by the Fick- Jacobs equation^*^'^'^ which is derived from the 
3D (or 2D) Smoluchowski equation after elimination of y and z coordinates by assuming 
equilibrium in the orthogonal directions. Reduction of the spatial dimensionality from two 
or three physical dimensions to an effective one dimensional system may involve not only 
the appearance of an entropic barrier, but also the effective diffusion coefficient. When 
\u {x)\ << 1, the effective diffusion coefficient 

where Dq = ksT/ri and a = 1/3 and 1/2 for two and three dimensions, respectively. The 
prime stands for the derivative with respect to the space variable x. 
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Consider the effective diffusion coefficient and the entropic barrier, the dynamics of a 
Brownian particle moving along the axis of the 3D (or 2D) tube can be described by^'^'^'^ 

dP{x,t) d , dP{x,t) , D{x)dA{x,t) ^^ dj{x,t) 

dt ~ dx^ ^""^ dx ^ ksT dx ^^^'^^J- dx ' ^ ^ 

where we define a free energy A{x^t) := E — TS = U{x) — F{t)x — T/cBln/i(x), here 
E = U{x) — F{t)x is the energy, S = kBlnh{x) is the entropy, h{x) is the dimensionless 
width 2uj{x)/L in two dimensions, and the dimensionless transverse cross section tt[u){x)/L]'^ 
of the tube in three dimensions. Here, we define a effective potential along the x coordinate, 

U^^f (x) = U{x)- TkB In h{x) . (9) 

t) is the probability current density. t) is the probability density for the particle 
at position x and at time t. It satisfies the normalization condition Jq P{x,t)dx = 1 and 
the periodicity condition P{x,t) = P{x + L,t). 

If F{t) changes very slowly with respect to t, namely, its period is longer than any other 
time scale of the system, there exists a quasisteady state. In this case, by following the 
method iii^"— , we can obtain the current 



= ^sTll-.M-^)] 



For the given external force (Eq. (4)), the average current^i^ is 

J=- rmt))dt = l[j{Fo)+j{-Fo)]. (11) 
T Jo 2 



III. RESULTS AND DISCUSSIONS 



Because the results from two and three dimensions are very similar, for the convenience of 
physical discussion, we now mainly investigate the current in three dimensions with /c^ = 1, 
1] = 1 and L = 27t. The results are shown in Figs. 2-5. 

Figure 2 shows the current contours on T — plane. There are two hills and two vales in 
the plot. The current is negative in area A and C, positive in area B and D. In our model, 
the potential and the tube shape in space and the external force in time are symmetric. 
It is obvious that the phase shift between the potential and the tube shape leads to a net 
current. However, the current is always zero at = 0, and 27r. When < < 0a and 
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(f)f) < (f) < 27r, the particle will undergo a hill and a vale on increasing temperature, namely, 
the current reverses its direction. When temperature is increased, the particle will undergo 
only a hill or a vale for 0a < < vr and n < (p < (pf,. The current is always zero at = vr. 
Because of the numerical trouble near = l.Ovr and low temperature, the numerical results 
at = l.Ovr and low temperature meet a little error. However, from the effective potential 
(Eq. (9)), It is easy to find that the effective potential is symmetric and the current is zero 
at = l.Ovr. So the line at = vr should be vertical. 

The current J as a function of temperature T is shown in Fig. 3 for different values of 
0. When T — )■ 0, the particle cannot pass the barrier, so the current J tends to zero. When 
T — )■ oo, the thermal noise is very large, so the ratchet effect disappears and the current J 
approaches zero, also. There exist optimized temperatures at which the current J takes its 
maximal or minimal values. When = 0.2vr (0 < < 0a) and 1.8vr(0;, < < 2vr), current 
reversals may occur on increasing temperature. The current is always positive and negative 
for = 0.7vr (0a < < vr) and 1.4vr (vr < < 0;,), respectively. It is easy to see from 
the Figs. 2 and 3 that the current on increasing temperature may reverse its direction for 
< < 0a and 0{, < < 2vr and is always positive at all temperatures for 0a < < vr or 
negative for vr < < 0;,. 

In order to illustrate the current (see Fig. 3) in detail, the effective potential Ueff{x) 
along the x coordinate is plotted in Fig. 4. From the previous works^"— , we can obtain 
that there is no current when the effective potential is symmetric. When the right of the 
potential in a period is steeper than the left, the current is positive. The current is negative 
when the left is steeper. In Fig 4 (a)(0 = 0.2vr), the left is steeper than the right (inducing 
a positive current) at low temperatures. However, when temperature is increased, the right 
may be steeper than the left (inducing a negative current) at high temperatures. For the 
case = 1.8vr (Fig. 4(d)), the current is negative at low temperatures and positive at high 
temperatures. Therefore, the current can change its direction on increasing temperature for 
these two cases (0 = 0.2vr and 1.8vr). In Fig. 4 (b) (0 = 0.7vr), the right is always steeper 
than the left, so the current is always positive. The left is always steeper than the right 
at = 1.4vr (Fig. 4 (c)), so the current is always negative. In these two cases, no current 
reversal occurs. 

Figure 5 shows the current J as a function of the phase shift 0. There is no net current 
at = 0, vr and 2vr for all temperatures. There are two peaks and two vales in the curve 
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on increasing the phase shift at T = 0.2 (low temperature). The current can change its 
direction even three times when the phase shift is changed. At high temperature (T = 0.5), 
only one peak and one vale appear and the current reverses its direction at = tt. 

IV. CONCLUDING REMARKS 

In present work, we study the transport of a Brownian particle moving along the axis of 
three-dimensional periodic tube in the presence of a temporally symmetric unbiased force 
and a symmetric potential. The movement equation can be described by the Fick- Jacobs 
equation which is derived from 3D (or 2D) Smoluchowski equation after elimination of y and 
z coordinates. Reduction of the coordinates may induce the appearance of an entropic barrier 
and the effective diffusion coefficient. The current is achieved by solving one-dimensional 
Fick- Jacobs equation. 

It is found that the phase shift between the tube shape and the potential may induce a net 
current even when the potential, the tube shape and the external force are symmetric. The 
symmetry of the effective potential can be broken by changing temperature or the phase shift. 
When < < 0a and (pb < (p < 27r, current reversal may occur on increasing temperature. 
However, the current is always positive for 0^ < < tt or negative for vr < < 0^. No 
current occurs at = 0, vr and 2tt for all temperatures. At low temperatures the current 
may reverse its direction three times on increasing the phase shift. At high temperatures 
the current changes its direction one time at = vr. 

Our tube model is to describe molecular motors, especially for kinesins and dyneins along 
microtubule. The potential describes the track of the motor. The external force F{t) depicts 
the stroke force due to ATP hydrolyzing. The shape of the tube describes the structure of 
the motor. The phase shift between the tube shape and the potential may induce a net 
current. Our results can explain current reversals in molecular motors. Several biological 
molecular motors, for instance kinesin and non-claret disjunctions, belonging to the same 
superfamily of motor proteins move towards opposite ends of the microtubule^-. This can be 
explained by current reversals 

The results we have presented have broad applications in many processes^^, such as 
molecular motors movement through the microtubule in the absence of any net macro- 
scopic forces^-, ion transport through ion channels^, motion of polymers subjected to rigid 



7 



The Journal of chemical physics 126, 204706 (2007) 



constrainta^^, drug release^^and polymer crystallization^i. In these systems, a directed- 
transport with entropic or energy barriers can be obtained in the absence of any net macro- 
scopic forces or in the presence of the unbiased forces. 
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VI. CAPTION LIST 

FIG. 1. The upper is schematic diagram of tube with periodicity L. The shape is 
described by the radius of the tube u{x) — asin(^ + 0) + 6. The bottom is the potential 
along the x coordinate with periodicity L, U{x) = Qsm(^). 

FIG. 2. Current contours on T — 0(7r) plane Bi Q — 0.5, a — b — a — 1/2 and 
F — 0.5. The current is negative in area A and C, positive in area B and D. The solid line 
denotes zero current. 0a (0.3817) and 0b (1.6035) are the intersectant point between the 
0-axis and the zero-current line near area A and D, respectively. 

FIG. 3. Current J vs temperature T for different values of at Q = 0.5, a—^,h — 
a = 1/2 and Fq = 0.5. 

FIG. 4. The effective potential Ueff{x) along the x coordinate for different values of T. 
From the bottom to the upper the temperature is increased, (a) = 0.27r; (b)0 — O.Itt; 
(c)0 = lAn; (d)0 = I.Stt. 

FIG. 5. Current J vs the phase shift for different temperature a,t Q — 0.5, a — 
h^]^^a = \l2 and Fq = 0.5. 
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